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ABSTRACT 

In this work we present Spitzer 3.6 and 4.5 /im secondary eclipse observations of five new cool (< 1200 
K) transiting gas giant planets: HAT-P-19b, WASP- 6 b, WASP-lOb, WASP-39b, and WASP-67b. We 
compare our measured eclipse depths to the predictions of a suite of atmosphere models and to eclipse 
depths for planets with previously published observations in order to constrain the temperature- and 
mass-dependent properties of gas giant planet atmospheres. We find that the dayside emission spectra 
of planets less massive than Jupiter require models with efficient circulation of energy to the night side 
and/or increased albedos, while those with masses greater than that of Jupiter are consistently best- 
matched by models with inefficient circulation and low albedos. At these relatively low temperatures 
we expect the atmospheric CH 4 /CO ratio to vary as a function of metallicity, and we therefore use 
our observations of these planets to constrain their atmospheric metallicities. We find that the most 
massive planets have dayside emission spectra that are best-matched by solar metallicity atmosphere 
models, but we are not able to place strong constraints on metallicities of the smaller planets in 
our sample. Interestingly, we find that the ratio of the 3.6 and 4.5 p m brightness temperatures for 
these cool transiting planets is independent of planet temperature, and instead exhibits a tentative 
correlation with planet mass. If this trend can be confirmed, it would suggest that the shape of 
these planets’ emission spectra depends primarily on their masses, consistent with the hypothesis that 
lower-mass planets are more likely to have metal-rich atmospheres. 

Subject headings: eclipses - planetary systems - techniques: photometric 


1. INTRODUCTION 

Observations of the thermal emission spectra of gas gi¬ 
ant planets provide an invaluable tool for probing their 
atmospheric compositions and pressure-temperature pro¬ 
files. To date nearly one hundred transiting planets have 
been observed in secondary eclipse with the Spitzer Space 
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Telescope , but the vast majority of these objects have 
been Jovian- mass p l anets with tem peratures between 


1500-2500 K (Madhusudhan et al. 2014). The atmo¬ 
spheres of smaller and cooler planets remain mostly un¬ 
charted territory, in part because it is extremely challeng¬ 
ing to detect their thermal emission at the near-infrared 
wavelengths accessible to most telescopes. Equilibrium 
chemistry models predict that between 1000-1200 K the 
dominant reservoir of atmospheric carbon shifts from CO 


atmospheres of brown dwarfs (Kirkpatrick 

2005 

). How- 

ever, these predictions are dependent not on 

ly on temper- 


ature but also on the underlying elemental abundances, 
and therefore are sensitive to changes in atmospheric 
metallicity. For cool planets where equilibrium chem¬ 
istry would normally predict a methane-dominated at¬ 
mosphere, a decrease in the relative hydrogen abundance 
will result in reduced CH 4 and enhanced CO abundances 


relative to the solar metallicity predictions (Moses et al. 
2013 Hu &; Seager|2014||Agundez et al.|20T4 |Venot et al. 


2014). 

This hypothesis was first proposed to explain the emis¬ 
sion spectrum of the Neptune-mass planet GJ 436b, 
which appears to have a CO-rich and methane-poor at- 



abundances for atmospheric metallicities up to 50 x so- 
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TABLE 1 

Target System Properties 


Target 

Stellar T e// (K) 

Mp (Mj up ) 

Period (days) 

T c a 

e 

Teg (K) b 

Ref. 

HAT-P-19b 

5007 ±66 

0.292 ± 0.018 

4.008778 ± 6 x 10"® 

5091.5342 ± 0.0003 

0.067 ± 0.042 

1010 

C 

WASP-6L 

5375 ± 65 

0.485 ± 0.027 

3.36100208 ± 31 x 10~ 8 

4425.02180 ±0.00011 

0.041 ± 0.019 

1184 

d 

WASP-lOb 

4735 ± 69 

3.14 ±0.27 

3.09272932 ± 32 x 10 -8 

4664.038090 ± 0.000048 

U.U4/3_q 0029 

980 

e 

WASP-39b 

5400 ± 150 

0.284 ±0.031 

4.0552965 ± 10 x 10 -6 

5342.96956 ± 0.00020 

0 (fixed) 

1116 

f 

WASP-67b 

5200 ± 100 

0.406 ± 0.035 

4.6144109 ± 27 x 10“ 7 

5824.37424 ± 0.00022 

0 (fixed) 

1004 

g 


a Measured time of transit center, BJDtjtc - 2,450,000. The reported times for WASP-6b, WASP-lOb, and WASP-67b are given in BJDtdb, and 
can be converted to the UTC time at the epoch of our observations by subtracting 66.2 s, XX s, and XX s, respectively. 
b Planet dayside equilibrium temperature calculated assuming zero albedo and uniform heat redistribution. 




|Gillon et al. 

(20091 

Husnoo et al. 

C2012|; 

Doyle et al.| (J2013||; 

Tregloan-Reed et al.1 (12015} 

ITorres et al. 

(2012) 

Barr os et al. 

2013); 

Knutson et al. (2014 

bT 

|Faedi et al. ( 

2011); 

Ricci et al. | (J 20151) 


Hellier et al. ( 2012 

;|Mancini et al. (2014| 


lar. Moses et al. (2013) explored atmosphere models 
spanning a broad range of metallicities and found that 
those with metallicities of 300 x solar or higher could pro¬ 
duce a good match to GJ 436b’s dayside emission spec¬ 
trum. Although a recent global re-a nalysis of all available 


Spitzer photometry for GJ 436b by Lanotte et ah ([20141 

Ste 


found a smaller eclipse depth at 3.6 /jm than Steven¬ 
son et ah, this study concludes that the revised emission 
spectrum is still best described by an atmosphere model 
with very little methane and significant CO. To date sec¬ 
ondary eclipse detections have been published for three 
additional gas giant planets with predicted equilibrium 
temperatures cooler than 1200 K, including WASP-8b 


(Cubillos et al. 2 013), WASP-80b ( Triaud et al. 
and HAT-F-20b(Deming et al. 2015), but none o 


2015), 

__ these 

planets appear to share GJ 436 b's unique atmospheric 
chemistry. Although HAT-P-12b was also observed with 
Spitzer, the secondary ec lipse was not detecte d in either 
the 3.6 or 4.5 /im bands ( Todorov et al.|[2013 ). 


Studies of hydrogen-dominated atmospheres in our 
own solar system have shown that as core mass fraction 
increases, so does the atmospheric metallicity. Uranus 
and Neptune have atmospheres with carbon to hydro¬ 
gen ratios approximately 70—100 times that of the solar 
value, while Jupiter’s atmospheric carbon to hydr ogen 


2004; Fletcher et al. 

2009) 

Karkoschka & Tomasko 2011| 

Sromovsky et ai.|201 

LI). A 

though we have relatively few 


trasolar gas giant planets, current mass and radius mea¬ 
surements suggest an inverse relationshi p between the 


pla net’s bulk metallicity and its mass (Miller & Fort¬ 
ney 12011 ). For sub-Neptune-sized planets, results from 


ney 

2011 ) 

the 

Keple 


tinue to increase with decreasing mass, with the caveat 
that small planets also appear to exhibit a greater diver¬ 
sity of densities at a fixed mass than their larger gas gi- 
ant c ounterparts (e.g. Marcy et al.||2014[ Weiss & Marc; 


201^Lqpez_^jFqrtri^ |p014[ Hadden &; Lithwick 2U 
Wolfgang fc Lopez||2015[ ). In this scenario we would not 

U/A Ul) Ou \ A A UG OAk U AT T) om. 4-^ ,, U--- 


f 


expect WASF-8b, WASP-80b, or HAT-P-20b to share 
Neptune-mass GJ 436b’s unique dayside chemistry, as 
these three planets all have masses greater than 0.5 M Jup . 

It currently remains an open question as to whether or 
not the metallicities of exoplanet atmospheres increase 
with decreasing mass, as models of planet formation in¬ 


dicate that small variations in protoplanetary disks can 
lead to a broad range of planet masses and compo sitions 
( Fortney et al.|20~T3 Helled fe Bodenheimer|20~T4 ). Mea¬ 
surements ot average density alone are not sufficient to 
uniquely determine atmospheric metallicity, as there is 
a degeneracy with the assumed com position and mass 
of the planet’s metal-rich core (e.g., Rogers fe Seager] 
20101 IFigueria et~ah |20091 INettelmann et al.| |2UlU| |Va-| 


lencia et al. 2013 Howe et al. 2014). Fortunately the 
combination of transmission spectroscopy and secondary 
eclipse observations provides unique leverage to directly 
constrai n the atmospheric meta llicities of transiting plan¬ 
ets (e.g., Kreidberg et al.|2014|. For planets with diffuse, 


high altitude cloud layers, secondary eclipse spectroscopy 
offers an added advantage as these cloud layers are less 
likely to be opti cally thick when viewed face-on than dur¬ 


ing transit (Fortney 2005). Similarly, planets with high 
mean molecular weight atmospheres may have absorp¬ 
tion features that are undetectable in transit with current 
telescopes. The Neptune-mass planet GJ 436b serves as 
a useful illustration of this concept, as it has a feature¬ 
less transmission spectrum but exhibits strong molecu¬ 
lar absorption features in its dayside emissio n spectrum 
( Stevenson et al.||2010 Knutson et al.|[2014c ). 


Observations ot transiting planets spanning a wide 
range of effective temperatures and masses can also re¬ 
veal correlations between the efficiency of atmospheric 
circulation and other properties of the system. Tidal evo¬ 
lution models predict that the tidal locking time scales 
for short period planets shoul d be significantly shorter 


than the ages of t hese systems (Guillot et al. 1996 Show¬ 


man et al. 2015). By measuring the day side temper¬ 


atures of a large, diverse sample of transiting planets, 
we can place constraints on the nature of the dynami¬ 
cal processes that transport heat between hemispheres, 
as well as their dependence on the incident stellar flux, 


vious studies by 

Cowan & Ago! (2011), Perez-Becker & 

Showman ( 

2013) 

, and SchwaA^& Cowan (2015) indicate 


that, on average, highly irradiated hot Jupiters appear 
to have inefficient day-night heat recirculation, whereas 
more weakly irradiated hot Jupiters seem to exhibit more 
efficient day-night circulation albeit with a greater degree 
of diversity. This trend is in good agreement with re¬ 
sults from three-dimensional general circulation model_s 
of hot Jupiters including radiative transfer (e.g., Show- 
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Fig. 1.— Raw Spitzer 3.6 and 4.5 pm photometry as a func¬ 
tion of time from the predicted center of eclipse for HAT-P-19b, 
WASP-6b, and WASP-39b. The relative flux from each star is 
shown binned in 30 s (grey filled circles) and 5 minute (black filled 
circles) intervals, and the best fit instrumental models are binned 
in 5 minute intervals and overplotted for comparison (solid lines). 
Intrapixel sensitivity variations cause distinct sawtooth patterns as 
a result of oscillations in the centroid position of the star with a 
period of approximately 45 minutes. Observations for HAT-P-19 
are plotted in chronological order, with visits from 2011-2012 one 
row up from the bottom and visits from 2014 in the bottom row. 


man et al. 2009, 2013, 2015; Heng 

2012; Rauscher & 

Menou 2012; Perna et al. 2012|) 

anc 

can be explained 

Perez-Becker & Show- 

with simple dynamical theories ( 


radiative transfer suggest that both the degree of devia¬ 
tion from radiative equilibrium on the planet’s day side 
and the overall magnitude of the infrared day-night flux 


on planetary rotation rate (ilShowman et a 

. 201 

5) and 

atmospheric metallicity ( 

Lewis et al. 2010 

. Al 

though 

there are multiple studies (Uowan & Agol 

2011 

Perez- 

Becker & Showman 2013 

Schwartz & Cowan 201 

5) that 


: explored ' 

for relatively hot planets, there are currently very few 
published observations constraining the efficiency of the 
day-night circulation on cooler planets as well as those 
with metal-rich atmospheres. 

In this work we present new Spitzer secondary eclipse 


Fig. 2. — Raw Spitzer 3.6 and 4.5 pm photometry as a function 
of time from the predicted center of eclipse for WASP-10b and 
WASP-67b; see Fig. [I] for more details. Duplicate observations for 
both planets are plotred in chronological order, with earlier visits 
in the upper row and later visits in the lower row. 


observations of the cool gas giant planets HAT-P-19b, 
WASP-6b, WASP-lOb, WASP-39b, and WASP-67b, 
which have predicted equilibrium temperatures between 
900—1200 K and masses of 0.3—3.1 Mj up . These observa¬ 
tions were obtained as part of a larger program designed 
to search for correlations between planet mass and atmo¬ 
spheric metallicity (GO 10054, PI Knutson), with a total 
sample size of fifteen planets. The selected planets have 
predicted equilibrium temperatures between 650-1150 K, 
and if they have hydrogen-rich atmospheres whose chem¬ 
istry is in local thermodynamic equilibrium we would ex¬ 
pect to see strong methane absorption features in their 
emission spectra. We list relevant characteristics for the 
five planets in this study in Table [I] In Section [2] we de¬ 
scribe the Spitzer data acquisition and reduction, while 
Sections [3] and [4] discuss the implications of these results 
for our understanding of the relationship between planet 
mass and atmospheric metallicity in the gas giant regime. 

2 . OBSERVATIONS AND DATA ANALYSIS 
2.1. Photometry and Instrumental Model 
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TABLE 2 

Spitzer Observation Details 


Target 

A (pm) 

UT Start Date 

Length (h) 

n img & 

tint (s) b 

i c 

L trvm 

n bin C 

rr> C 

1 pos 

c 

1 npzx 

Fixed d 

r phot° 

Bkd (%) e 

HAT-P-19b 

3.6 

2012-03-12 

7.9 

14144 

2.0 

0.5 

128 

2.5 

none 

yes 

2.2 

1.88 


3.6 

2014-10-02 

6.1 

10816 

2.0 

0.5 

16 

1.5 

2.5 

no 

1.4 

0.91 


4.5 

2011-09-10 

7.9 

2167 

12.0 

0.5 

8 

2.5 

2.5 

no 

2.5 

0.82 


4.5 

2014-09-28 

6.1 

10816 

2.0 

0.0 

32 

1.5 

2.5 

no 

1.4 

0.65 

WASP-6b 

3.6 

2009-12-30 

7.6 

2096 

12.0 

1.5 

8 

1.5 

none 

yes 

2.3 

-0.38 


4.5 

2010-01-23 

7.6 

2096 

12.0 

0.5 

5 

1.5 

none 

yes 

2.1 

1.47 

WASP-lOb 

3.6 

2011-01-19 

7.8 

3976 

6.0 

0.5 

64 

1.5 

none 

yes 

2.1 

-0.17 


3.6 

2014-09-10 

4.9 

8704 

2.0 

0.5 

64 

1.5 

1.9 

no 

2.2 

1.47 


4.5 

2010-09-17 

7.8 

3976 

6.0 

0.5 

64 

1.5 

none 

yes 

2.7 

-0.18 


4.5 

2014-09-22 

4.9 

8704 

2.0 

0.5 

4 

1.5 

none 

yes 

2.3 

0.39 

WASP-39b 

3.6 

2012-04-04 

7.0 

12480 

2.0 

2.0 

32 

2.5 

3.0 

no 

1.8 

1.29 


4.5 

2012-04-08 

7.0 

12480 

2.0 

1.0 

64 

1.5 

none 

yes 

2.4 

1.20 

WASP-67b 

3.6 

2014-06-28 

6.0 

10752 

2.0 

0.5 

8 

1.5 

none 

yes 

2.2 

1.49 


3.6 

2014-07-21 

6.0 

10752 

2.0 

1.0 

8 

1.5 

none 

yes 

2.3 

2.10 


4.5 

2014-07-07 

6.0 

10752 

2.0 

0.5 

16 

2.0 

2.1 

no 

1.3 

0.63 


4.5 

2014-07-30 

6.0 

10752 

2.0 

0.5 

64 

2.0 

none 

yes 

3.0 

0.97 


a Total number of images, 
b Image integration time. 

c ttrim is the amount of time in hours trimmed from the start of each time series, ribin is the bin size used in the photometric fits, r pos is the 
radius of the aperture used to determine the position of the star on the array, r np i x is the radius of the aperture used to calculate the noise pixel 
parameter, and r p hot is the radius of the photometric aperture (we provide the median aperture radius over the observation for the time-varying 
aperture case). All radii are given in pixels. 

b Denotes whether the photometry was obtained using a fixed or time-varying aperture. 
e Sky background contribution to the total flux measured in the selected aperture. 


These observations were obtained in the 3.6 and 4.5 fan 
bandpasses using the Infra-Red Array Camera (IRAC) 
on the Spitzer Space Telescopes. HAT-P-19 was ob¬ 
served in full array mode with no peak-up for the 2011 

4.5 /im eclipse, and in subarray mode with peak-up point¬ 
ing (this pointing adjustment corrects the initial tele¬ 
scope pointing in order to place the star near the cen¬ 
ter of the pixel) for all subsequent visits. The 2012 

3.6 fan eclipse was observed with a single peak-up at 
the start of the observation while the 2014 observations 
used the current standard approach, which includes both 
an initial peak-up of the position of the star on the ar¬ 
ray and a second peak-up thirty minutes into the obser¬ 
vations in order to correct the initial pointing drift of 
the teleOscope. WASP-6 was observed in full array mode 
with no peak-up pointing, as this mode was not available 
when these data were taken. WASP-10 was observed in 
full array mode with no peak-up pointing in the initial 
2010 — 2011 observations, while the 2014 observations 
were obtained in subarray mode with the standard dual 
peak-up pointing adjustments. WASP-39 was observed 
in subarray mode with no peak-up pointing, and WASP- 
67b was observed in subarray mode with the standard 
dual peak-up pointing method. Additional observation 
details can be found in Table [2] 

For each data set, we calculate BJDutc mid-exposure 
times using information in the BCD files provided by 
the Spitzer pipeline. We estimate and subtract the sky 
background, calculate the flux-weighted centroid position 
of the star on the array, and calculate the total flux in a 
circular aperture for each individual image as described 
in previous studies (IKnutson et al. 2012 Lewis et al. 


2013 Todorov et al. 2013 hammer et al. 2011). We show 


the resulting light curves in b ig. [I] and [5] We consider 
both fixed and time varying photometric aperture sizes in 
our fits. For the fixed apertures we consider radii between 
2.0 — 5.0 pixels, where we step in increments of 0.1 pixels 
between 2.0 — 3.0 pixels and in 0.5 pixel increments for 


larger radii. The time varying apertures were calcul ated 
based on a scaling of the noise pixel parameter (Mighell 


2005), which is proportional to the square of the full 
width half max of the stellar point spread function, and 
described by Eq. |T] below: 


(E4) : 

n 


( 1 ) 


where I n is the measured intensity of the n th pixel. 
We iteratively re-scale the noise pixel aperture radii as 
r = a\fj3 + C, where a is a scaling factor with a value 
of [0.6, 0.7, 0.8, 0.85, 0.9, 0.95, 1.0, 1.05, 1.1, 1.15, 1.2], 
and C is a constant with a value between —0.8 and +0.4 
stepping in 0.1 pixel increments. 

We account for variations in intrapixel sensitiv ity using 
a pixel -leve l decorrelation (PLD) method, which Deming 
et al. ([2015]) found to produce results that are superior 
to those" oF a simple polynomial decorrelation or pixel 
mapping method. Details of this PLD technique are de¬ 
scribed in Deming et al., but the basic methodology is 
as follows. First, we obtain the raw flux values for a 
3x3 grid of pixels centered on the position of the star, 
and then normalize these individual pixel values by di¬ 
viding by the total flux in each 3x3 postage stamp. 
This normalization effectively removes any astrophysical 
signals, and we expect that the remaining variations in 
the fluxes of individual pixels as a function of time are 
primarily due to changes in the position of the star on 
the array. Equation[2]shows the resulting forward model, 
which we fit simultaneously with the eclipse model: 




S'model {l) — 


( 2 ) 


e m 

i 

where F mo d e i is the predicted stellar flux in an individual 
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TABLE 3 

Best Fit Eclipse Parameters 


Target 

Band (/xm) 

Fp/F* {%) 

F P /F,, avg (%r 

Tbright (K) a 

T s h 

O-C (d) c 

e cos uj 

HAT-P-19b 

3.6 

0.071 ± 0.020 

0.062% ±0.014% 

1090±« 

5999.5293 ± 0.0074 

0.0062 ± 0.0075 

0.0015 ±0.0021 


3.6 

0.053 ±0.021 


6933.5697 ± 0.0071 

0.0013 ±0.0076 



4.5 

0.066 ±0.019 

0.062% ± 0.016% 

9141“ 

5815.1239 d 




4.5 

0.053 ±0.028 


6929.5614 d 



WASP-6b 

3.6 

0.094 ± 0.019 


12351™ 

5196.3606 ±0.0106 

-0.0110 ±0.0106 

0.0001 ± 0.0020 


4.5 

0.115 ±0.022 


1118+®® 

5219.9007 ± 0.0044 

0.0021 ± 0.0044 


WASP-lOb 

3.6 

0.103 ±0.015 

0.100% ±0.011% 

H52l H 

5580.9226 ± 0.0025 

-0.1095 ±0.0025 

-0.0552 ±0.0007 


3.6 

0.097 ± 0.017 


6910.7917 ±0.0027 

-0.1141 ± 0.0027 



4.5 

0.126 ±0.025 

0.146% ±0.016% 

1086139 

5457.2249 ± 0.0040 

-0.0980 ± 0.0040 



4.5 

0.158 ±0.020 


6923.1693 ±0.0020 

-0.1077 ± 0.0020 


WASP-39b 

3.6 

0.088 ±0.015 


1213l|® 

6022.2356 ± 0.0066 

0.0033 ± 0.0068 

0.0007 ± 0.0017 


4.5 

0.096 ±0.018 


10551“ 

6026.2878 ± 0.0058 

0.0005 ± 0.0060 


WASP-67b 

3.6 

0.014 ± 0.018 

0.022% ± 0.013% 

< 1046 s 

6860.3124 d 




3.6 

0.032 ±0.020 



6837.2416 d 




4.5 

0.088 ± 0.024 

0.080% ± 0.018% 

10151“ 

6846.4689 ± 0.0035 

0.0026 ±0.0041 

0.0007 ± 0.0012 


4.5 

0.070 ± 0.026 


6869.5390 ± 0.0064 

0.0007 ± 0.0068 



a We report the error-weighted mean eclipse depths at 3.6 and 4.5 for planets with multiple observations in each bandpass. We note that the 
averaged 3.6 /im eclipse depth for WASP-67b is formally a non-detection with a significance of 1.7a. Brightness temperatures are calculated using a 
PHOENIX stellar model interpolated to match the published temperature and surface gravity for each star. 

k BJDutc - 2,450,000. To convert from UTC to TDB time standards, simply add 66.2 s to eclipse times obtained between 2009 — 2012, and 67.2 s to 
eclipse times from 2014. 

c Observed minus calculated eclipse times, where we have accounted for the uncertainties in both the measured and predicted eclipse times as well 
as the added delay from the light travel time in the system. This delay is 61 s, 54 s, 49 s, 63 s, and 66 s for HAT-P-19b, WASP-6b, WASP-lOb, 
WASP-39b, and WASP-67b, respectively. 

^ We allow the eclipse times in this bandpass to vary as free parameters in our fit, but we use the error-weighted mean orbital phase and corresponding 
uncertainty from the other bandpass to place a prior constraint on their values. This red uces the uncertainty in the measured eclipse depths for visits 
where the eclipse depth is not detected at a statistically significant level (> 3a; see |2.2| for more details). 
e 2a upper limit based on the error-weighted average of the two 3.6 /im eclipse measurements. 


image, Fi is the measured flux in the i th individual pixel, 
and Wi is the weight associated with that pixel. We leave 
these weights as free parameters in our fit, and solve for 
the values that best match our observed light curves. 

We consider versions of the photometry with varying 
bin sizes, aperture sizes, and with up to two hours of data 
trimmed from the start of the time series in our fits. By 
fitting the binned light curves we are able to identify so¬ 
lutions with significantly less red noise on longer time 
scales and slightly higher scatter on the shortest time 
scales. The eclipses of the planets in this study typi¬ 
cally have durations of several hours, and we are there¬ 
fore most sensitive to noise on these time scales in our 
fits. For each fit we determine the best-fit model solution 
using the binned photometry, and then apply our best-fit 
model to the unbinned light curves in order to generate 
the plots shown here. We calculate our best-fit eclipse 
model and corresponding measurement uncertainties us¬ 
ing the version of the photometry that minimizes the 
amount of time-correlated (“red”) noise in the residuals 
after the best-fit model has been removed. As discussed 


in Derning et al. (2015), we evaluate the noise properties 
of a given set of residuals by calculating the root mean 
variance of the residuals as a function of bin size, where 
we step through bin sizes in intervals of powers of two. If 
the noise is Gaussian and there is no correlation between 
adjacent data points, we would expect the RMS of the 
binned data to decrease as 1/ 1 Jn where n is the number 
of points in each bin. We evaluate the relative amount 
of red noise in each version of the photometry by tak¬ 
ing the difference between a model with the ideal 1 /y/n 
scaling and the observed RMS at each bin size, squaring 


the difference, and summing over all bins in order. By 
minimizing this least-squares metric, we can select the 
version of the photometry which has the least amount of 
red noise across all time scales (see Fig. [3] and |4j) . 

When considering the optimal choice of aperture and 
bin size, we first identify the version of the photome¬ 
try with the lowest overall RMS variance in the best-fit 
residuals and discard all of the alternative versions of 
the photometry with a RMS variance that is a factor of 
1.2 greater than this value. This ensures that we avoid 
solutions that minimize the red noise component by sig¬ 
nificantly increasing the amount of white noise (these are 
usually the largest apertures, where the sky background 
contributes additional noise). We then pick the bin size 
and photometric aperture that has the smallest amount 
of red noise as measured by our least squares metric. 

We find that there is an exponential ramp visible at 
the start of some visits, which is a well-known feature 
of the IRAC 3.6 and 4.5 /jm arrays and has a shape 
that varies depending on the illumination history of the 
array immediately prior to each obse rvatio n as well as 
the brightnes s of the target star (e.g., Lewis et al.||2013~ 


Zellern et al. 2014). We remove this ramp by trimming 
data from the start of each light curve. After identifying 
the optimal bin size and aperture choice, we determine 
the correct amount of data to trim by examining the best 
light curve after detector effects have been removed and 
trimming until no ramp is visible at the start of the ob¬ 
servations. We then revisit our choice of bin size and 
aperture with the new trim duration in order to make 
sure that our previous choice is still the best one with 
the new trim interval. As shown in Fig. [3] and [4j this 
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Fig. 3.— Standard deviation of the residuals for HAT-P-19, 
WASP-6, and WASP-39 after the best-fit eclipse and instrumental 
noise models have been subtracted as a function of bin size (black 
lines). We also show the expected 1 /y/n scaling for Gaussian mea¬ 
surement errors and no correlation between adjacent photometric 
points (white noise) as red dashed lines, where we have normal¬ 
ized both curves to match the standard deviation of the unbinned 
residuals. Observations for HAT-P-19 are plotted in chronological 
order, with visits from 2011-2012 one row up from the bottom and 
visits from 2014 in the bottom row. 


results in light curves with no detectible red noise com¬ 
ponent on the time scales spanned by our observations. 
We also calculate a Lomb-Scargle periodogram for the 
best-fit residuals from our nominal light curve for each 
visit and find that for periods between five minutes and 
five hours there are no peaks in the power spectrum with 
false alarm probabilities lower than 10%. As a final step, 
we check to make sure that we obtain eclipse depths that 
are consistent at the ler level across different versions of 
the photometry that meet our 20% excess noise threshold 
and have minimal red noise components. 


2.2. Eclipse Model and Uncertainty Estimates 


We use the routines described in Mandel & Agol (2002) 
in order to generate our eclipse light curves, where we set 
the planet-star radius ratio, orbital inclination, and the 
ratio of the orbital semi-major axis to the stellar radius 
equal to their best-fit values from previous studies (see 
Table [l] for a complete list of references for each planet). 


Fig. 4.— Standard deviation of the residuals for WASP-10 and 
WASP-67 after the best-fit eclipse and instrumental noise models 
have been subtracted as a function of bin size. Observations for 
both planets are plotted in chronological order, with earlier visits 
in the upper row and later visits in the lower row. See Fig. [3] for 
more details. 


We assume a circular orbit for the purpose of calculating 
the eclipse shape and duration, as there is no evidence for 
an orbital eccentricity large enough to appreciably alter 
the eclipse duration (e sin uj > 0.1, where e is the orbital 
eccentricity and w is the longitude of periastron) in any 
of these systems. We allow the individual eclipse depths 
and center of eclipse times to vary as free parameters 
in our fits. For the second 4.5 /im eclipse observations 
of HAT-P-19b and both 3.6 /im observations of WASP- 
67b the eclipse is detected with marginal significance and 
we therefore place a Gaussian prior on the phase of the 
secondary eclipse in order to reduce the corresponding 
uncertainties on the eclipse depth. We do not fix the 
center of eclipse time in these fits, but instead implement 
this prior as a penalty in % 2 proportional to the deviation 
from the error-weighted mean center-of-eclipse phase and 
corresponding uncertainty from the two 3.6 (HAT-P-19b) 
or 4.5 (WASP-67b) /im eclipses. 

We find that the use of this prior changes the best-fit 
eclipse depth for the second 4.5 /im observation of HAT- 
P-19b by less than O.ler, while reducing the measurement 
error on this depth by 40% as compared to a fit with no 
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Fig. 5.— Normalized Spitzer 3.6 and 4.5 pm light curves for 
HAT-P-19, WASP-6, and WASP-39 as a function of time from the 
predicted center of eclipse, where we have divided out the best-fit 
instrumental model shown in Fig. The normalized flux is binned 
in 10 minute intervals, and best nt eclipse model light curves are 
over plotted for comparison (solid lines). Observations for HAT- 
P-19 are plotted in chronological order, with visits from 2011-2012 
one row up from the bottom and visits from 2014 in the bottom 
row. 


constraint on the center of eclipse phase. For consis¬ 
tency, we apply the same prior to the first 4.5 pm eclipse 
observation of HAT-P-19b, which is detected with a sig¬ 
nificance of 3.2d in fits with the standard uniform prior 
on the eclipse phase. We find that in this fit the use of 
the prior has no effect on the best-fit eclipse depth and 
reduces the measurement error in this parameter by 5%, 
as expected for a case where the data themselves provide 
good constraints on the eclipse time. For WASP-67b our 
3.6 /im fits with a uniform prior prefer an eclipse solu¬ 
tion at a much earlier phase with a best-fit eclipse depth 
that is consistent with zero at the lcr level. We there¬ 
fore conclude that this is likely a local minimum in y 2 
caused by a trade-off with the instrumental noise model, 
and place a prior on the eclipse phase based on the mea¬ 
sured 4.5 pm eclipse phase. This allows us to obtain new 
constraints on the eclipse depth at 3.6 /(m for a range of 
phases consistent with the measured 4.5 /im eclipses. 

We determine the best fit model for each individual 


Fig. 6.— Normalized Spitzer 3.6 and 4.5 pm light curves for 
WASP-10 and WASP-67 as a function of time from the predicted 
center of eclipse, where we have divided out the best-fit instru¬ 
mental model shown in Fig. [2] Observations for both planets are 
plotted in chronological order, with earlier visits in the upper row 
and later visits in the lower row. See Fig. [5] for more details. 

eclipse using a Levenberg-Marquardt minimization rou¬ 
tine, where each fit has thirteen free parameters. In addi¬ 
tion to the nine pixel weight parameters and two eclipse 
parameters, we also fit for a linear function of time in 
order to account for long-term instrumental and stellar 
trends. Figures]!] and [2] show the raw photometry, while 
Figures [5] and [BTshow the corresponding normalized pho¬ 
tometry after dividing out the best-fit PLD model and 
linear function of time. We estimate the corresponding 
uncertainties on our model parameters using a Markov 
chain Monte Carlo (MCMC) analysis with 10 6 steps per 
chain, where we define the lcr uncertainties as the sym¬ 
metric interval around the median parameter value con¬ 
taining 68% of the total probability. We assume uni¬ 
form priors on all parameters with the exception of the 
3.6 /cm eclipse phase of WASP-67b and the 4.5 /cm eclipse 
phase for HAT-P-19b as noted in the previous paragraph. 

We initialize our chains using the best-fit solution from 
the Levenberg-Marquardt minimization in order to min¬ 
imize the burn-in time. We then check to see where the 
X 2 value of the chain first drops below the median value 
over the entire chain, and trim all points prior to this step 
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Fig. 7.— Posterior probability distributions and covariances with 1, 2, and 3c contour levels for the secondary eclipse depth, observed 
minus calculated center of eclipse time, and linear function of time from from our MCMC analysis of the WASP-6b 3.6 pm observation. 
We calculate the O-C eclipse time assuming that the planet has a circular orbit and using the transit ephemeris from |Tregloan-Reed et ah] 
(2015). We find no correlation between the measured eclipse depth and the other two parameters, but there does appear to be a correlation 
between the center of eclipse time and the linear trend slope, which is likely responsible for the larger than usual measurement errors on 
our reported center of eclipse time. 


in order to remove any residual burn-in phase. We also 
examine our resulting posterior probability distributions 
in order to make sure that our choice of prior bounds is 
broad enough that it does not affect our final results, and 
re-run our chain with a broader prior range if required. 
Although we find that individual pixel weights are of¬ 
ten degenerate with each other, our probability distribu¬ 
tions for the best-fit eclipse depths and center of eclipse 
times are Gaussian and uncorrelated with the other fit 
parameters, with the exception of the 3.6 /im eclipse of 
WASP-6b which we discuss in more detail below. The 
correlation between pixel weights could be mitigated by 
a tra nsformation to an independ ent set of basis vectors 


(e.g., Morello et al. 2014 20151, but as these correla¬ 


tions do not affect our final best-fit eclipse parameters 
and their corresponding uncertainties we do not pursue 
this option here. 

We find that our measured center of eclipse time for 
WASP-6b’s 3.6 /un eclipse has an uncertainty that is sig¬ 
nificantly larger (15 minutes) than is typical for our other 
observations with statistically significant eclipse detec¬ 
tions. We examined the results from our MCMC analysis 
and identified a correlation between the slope of the lin¬ 
ear trend across the 3.6 micron eclipse and the measured 
center of eclipse time in this visit (see Fig. [7]). As dis¬ 
cussed earlier, we found it necessary to trim the first 1.5 
hours of data in order to remove an unusually long-lived 
ramp at the start of these observations. This trimming 
means that we have a relatively short baseline prior to 


the start of the eclipse, and it is therefore not surpris¬ 
ing that there might be a correlation between the eclipse 
time and the linear trend in our fits. Our best-fit eclipse 
depth for this visit does not appear to be correlated with 
either of these parameters, and we find no evidence for a 
similar correlation between eclipse time and linear trend 
slope in other visits with short pre-eclipse baselines. 

3. DISCUSSION 
3.1. Orbital Eccentricity 

WASP-lOb is the only planet in our study that was 
previously known to have an eccentric orbit, with e = 
0.0473 ± 0.0032 based on a total of sixteen radial ve- 
locity measurements from Keck HIRES a nd SOPHIE 
( Christian et al.|2009 Knutson et al.|2014b| ). Our error- 
weighted average secondary eclipse phase of 0.4649 ± 
0.004 confirms this planet’s non-zero orbital eccentricity, 
and we combine ou r measurement with pub lished radial 
velocity data from Knutson et al. (2014b) in order to 
obtain improved constraints on the orbital eccentricity 
e and longitude of periastron w (see Fig. |8j. We also 
add two new Keck HIRES measurements that extend 
the baseline of these data by approximately a year and 
a half beyond the last measurement reported in our pre¬ 
vious paper. These new measurements clearly establish 
the presence of curvature in the observed radial velocity 
trend, and we have adapted our fits accordingly. 

We give the updated orbital parameters from a joint fit 
to the radial velocity measurements and using the mea- 
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Fig. 8 . — The top panel shows radial velocity measur e ments 
(fill ed black circles) for WA SP-10 from IChristian et al.| (2009) 
and |Knutson et al.| (|2014b|) , with best-fit model over plotted in 
blue. The residuals after subtracting the orbit of the transiting hot 
Jupiter WASP-10b are shown in the middle panel, with best-fit ra¬ 
dial velocit y trend overplott ed in re d. Although a linear trend was 
reported in [Knutson et al.| (|2014b|), the addition of two new Keck 
HIRES radial velocity measurements required us to fit a curved 
trend in this analysis. The best-fit orbital solution for WASP-lOb 
is shown in blue in the lower panel, with phased radial velocity 
measurements after re moval of the best-fit tre nd over plotted as 
filled black circles. See iKnutson et al.l (|2014b|) for additional de¬ 
tails. 


sured trans it ephemeris and seco ndary eclipse times in 
Table[4] (see|Knutson et al.j2014b for additional informa¬ 
tion on the fit 


itting method). Our new orbital eccentricity 
of e = 0.0589lg oQ 2 g differs from our previous value by 
4a; this change can be attributed to the inclusion of our 
measured secondary eclipse times in the new fits as well 
as the switch to a curved radial velocity trend model. We 
note that a previous study reported a detection of the 
secondary eclipse for this planet in the Ks band at an 
orbital phase of 0.4972 ± 0.0005 (Cruz et al. 2015); this 
offset differs significantly from that of our Spitzer sec¬ 
ondary eclipse observations and is also inconsistent with 
the current radial velocity measurements for this system. 
The authors of this study note that their reported error 
is likely an underestimate as it does not account for ad¬ 
ditional uncertainties contributed by the removal o f sys- 
tematic noise sources from their light curve. See §3.2. 1| 
for additional discussion of this observation. 

Published radial velocity observations for the remain¬ 
ing four planets in this study suggest that their orbits 
are likely circular, in good agreement with our measured 
secondary eclipse times. WASP-6b has the most exten¬ 
sive radial velocity data set, with a total of 79 measure¬ 
ments of which 38 were obtained in a single night in order 
to measure the Rossiter-McLaughli n effect as the p lanet 
passed in front of its host star (Gillon et al. 2009). Al¬ 
though the initial discovery paper reported a non-zero 


TABLE 4 

Updated Orbital Solution for WASP-10 from 
Radial Velocity Fit 


Parameter 

Value 

Units 

p b 

3.0927278 ±2.8e - 06 

days 

T 

-^conj ,b 

2454664.03809 ±4.7e - 05 

BJD tdb 

K b 

556.7 tit 

m s —1 

ecos uj 

-0.054 ±0.001 


e sin a; 

0.027 tlZ 1 


eb 

0.0608 1" 


UJ b 

153.6 tlr 

degrees 

71 

15.3 ±2.7 

m s —1 

72 

-64 ±41 

m s —1 

7 

-0.037 ±0.013 

m s —1 day —1 

7 

-5.94e-05 

m s -1 day -2 

jitter 

6.3 

m s —1 


These fits used the transit ephemeris from ^Barros^e^jilJ 
j2013fr and the secondary eclipse times reported in this pa- 
per as priors on the fit in order to obtain improved con¬ 
straints on the orbital solution. Radial velocity data taken 
from [Christian et al. | ( |2009| ) and |Knutson et al. | ( |2014b| ) , with 
the addition of two new unpublished RV measurements from 
Keck/HIRES. We allow the zero points of each radial veloc¬ 
ity data set (71 and 72 to vary independently in our fits, 
and describe the radial velocity trend with a linear ( 7 ) and 
quadratic ( 7 ) function of time. The reference epoch for 7 , 7,7 
is 2455615.0 


orbital eccentricity of 0.054 ±0.017, th is find ing was sub¬ 
sequently disputed by Husnoo et al-1 (|2012). These au¬ 
thors reanalyzed the same data while excluding the sin¬ 
gle night corresponding to the Rossiter measurement and 
found a best-fit orbital eccentricity of 0.04 ± 0.02. For 
HAT-P-19b the radial velocity data constrain the orbital 
eccentricity to a value of 0.07 ± 0.04, and for WASP-39b 
the orbi tal eccent ricity was fixed to zero in the published 
fits (Faedi et al. 2011). Spitzer secondary eclipse obser- 
vations provide a significantly tighter constraint on the 
value of ecosta, where e is the orbital eccentricity and ca 
is the longitude of periastron. These eclipse times can be 
used to further refine the global o rbital eccentricity when 
fitting radial velocity data (e.g., |Knutson et al.||2014b|). 
We find that the observed secohdary eclipse times tor 
these four planets are consistent with circular orbits, and 
list the corresponding ecosca values in Table [3] 

3.2. Comparison to Atmosphere Models 

We next compare the measured Spitzer secondary 
eclipse depths for each planet to the predictions of stan¬ 
dard one dimensional atmospher e models. We c o nsider 
two classes of mod e ls, ba sed on |Burrows et al.| (2008) 
and Fortney et al. (2008). Both models assume that 
the chemistry of these atmospheres is in local thermo¬ 
dynamic equilibrium and parameterize the unknown re¬ 
circulation of energy to the night side. 


3.2.1. Constraints on Atmospheric Metallicity and 
Day-Night Circulation 


We use the Fortney et al. (2008) models to explore 
a range of atmospheric metailicities from 1 — 1000 x so¬ 
lar. For planets at these temperatures, an increase in the 
atmospheric metallicity results in a decrease in the rel¬ 
ative amount of methane as compared to CO and CO 2 . 
Because methane absorbs in the 3.6 gm Spitzer band- 
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Fig. 9.— Predicted planet-star flux ratio as a function of wave¬ 
length for lx solar (purple) and lOOOx solar (light blue) atmo¬ 
spheric metallicities and e fficient recirculation of energy from the 
day side to the night side (|Fortney et al.|2008). We also show the 
corresponding lx solar (red) and lUUUx solar (orange) models for 
inefficient recirculation. Both models assume a dayside albedo of 
zero, where a non-zero albedo would have the same effect as increas¬ 
ing the efficiency of the day-night circulation. Models are shown 
as solid lines, with the corresponding band-integrated model pre¬ 
dictions overplotted as filled squares. Measured Spitzer secondary 
eclipse depths with their associated 1 a uncertainties are shown as 
black filled circles, with the exception of the 3.6 pm eclipse of 
WASP-67b where we plot the 2cr upper limit. In addition to the 
three planets described in this paper, we also compare th ese same 
models to p ublished secon dary eclipse depths fo r GJ 436b d Lanotte 
et al.||2014|), HAT-P-20b (Doming et al. 2015|, WA SP-8b dOubTF 
ios et al.|2(Jl3l iDeming et ai.|2U15|), and WAbP-80b flTriauaet al. 
2U15|). We list masses and predicted equilibrium temperatures tor 
each planet at the time of secondary eclipse calculated assuming 
zero albedo and efficient redistribution of energy to the night side. 


pass while CO and CO 2 absorb in the 4.5 /im band¬ 
pass, metal-rich models will have correspondingly strong 
emission at 3.6 /im and weak emission at 4.5 fi m as com- 



Fig. 10.— Predicted planet-star flux ratio as a function of wave¬ 
length for lx solar and lOOOx solar atmospheric metallicities and 
either efficient or inefficient recirculation of energy from the day 
side to the night side (see Fig. [9]for more information). Measured 
Spitzer secondary eclipse depths with their associated la uncer¬ 
tainties are shown as black filled circles, with the exception of the 
4.5 pm eclipse of GJ 436bb where we plot the 2a upper limit. 


pared to their solar metallicity counterparts. At 1 x solar, 
the pressure-temperature profiles in the models presented 
here are in radiative equilibrium and are self-consistent 
between the predicted equilibrium chemistry mixing ra¬ 
tios, the opacity of these molecules, and the the absorp¬ 
tion and emission of flux. At lOOOx solar, the models 
use the radiative-equilibrium pressure-temperature pro¬ 
files from a 50 x solar model where we have recalcu¬ 
lated the chemical equilibrium and spectra at lOOOx so¬ 
lar metallicity. Although not fully self-consistent, these 
models are a reasonable approximation for the appear¬ 
ance of a very metal-enhanced atmosphere. We test the 
validity of these models in the high metallicity regime 
by comparing them t o the measur ed sec ondary eclipse 
depths for GJ 436b (Lanotte et al. 2014). This planet 
appears to be best-matched by a model with lOOOx so¬ 
lar metallicit y in good agreement with previous studies 
(e.g., Moses et al. 2013), but we are not able to fully 
match the steep slope across the 3.6 and 4.5 pm bands. 
We expect that models including disequillibrum chemical 
process and/or non-solar carbon to oxygen ratios (e.g., 







































11 


Visscher fe Moses|2011 iMoses et al.|2011 ) could provide 
an improved match to these observations, and leave this 
as an exercise fo r fut ure modeling studies. 

Figure [9] and [lO] compare our measured secondary 
eclipse depths for the five planets in this study to lx 
and lOOOx solar metallicity models with either efficient 
or inefficient redistribution of energy to the night side. 
The relative efficiency of recirculation is implemented by 
varying the amount of flux incident at the top of the one¬ 
dimensional atmospheric column, and is therefore equiv¬ 
alent to a change in the planet’s dayside albedo. For 
the inefficient case, the incident flux is set equal to the 
dayside average value, while for the efficient case the inci¬ 
dent flux is calculated as the average over the entire sur¬ 
face. For four of the five new planets, our observations 
are best-matched by models with efficient recirculation 
of energy to the night side but we are unable to distin¬ 
guish between low and high metallicity atmosphere mod¬ 
els with high statistical significance. WASP-lOb appears 
to be the exception to this trend, and is best matched by 
models with inefficient recirculation and a relatively low 
atmospheric metallicity. It is also by far the most massive 
planet in our sample, a fact that we discuss in more detail 
below. Although we do not show it in these plots, there is 
a published K s (2.1 /im) secondary eclipse measurement 
available for this planet with a value of 0.137% ± 0.016% 


(Cruz et al. 


20151. The measured planet-star flux ratio 


in this bandpass is 3 — 10 x larger than the predictions 
of our models, and the reported center of eclipse phase 
differs from our measurement by 2.4 hours (54<r), compa¬ 
rable to the eclipse duration. It is therefore unlikely that 
this ground-based measurement corresponds to the same 
eclipse signal detected in our Spitzer photometry, and in 
fact the authors note that their reported uncertainties 
are likely underestimated due to the unknown contribu¬ 
tion of systematic noise sources to their error budget. 

Given that we obtain relatively tight constraints on 
the atmospheric properties of WASP-lOb, it is interest¬ 
ing to ask whether or not it is representative of other 
Jovian planets in this temperature range. A search of 
the literature indicates that there are two other planets 
(HAT-P-20b and WASP-8b) with measured Spitzer sec¬ 
ondary eclipse depths at 3.6 and 4.5 pm, masses greater 
than that of Jupiter, and predicted equilibrium temper¬ 
atures less than 1200 K. We compare the observations 
for these two planets to the same set of models as before, 
and plot the results in Figure [9] and [lO] For HAT-P- 
20 b the data are best matched oy the solar metallicity 
model with inefficient circulation of energy_to t he n ight 


(J2015 

and 

Triaud et al. (12015). 

et al. 

2013 

Dennng et al. 2fll5 


with inefficient energy recirculation and solar metallic¬ 
ity, but this model still under-predicts the measured flux 
in the 3.6 p m band. The unusual nature of WASP-8b’s 
daysi de emi ssion spectrum is discussed in detail in |Cu-| 
billos et al. (2013), and may be related to its eccentric 
orbit (e = 0.30). The only other eccentric planets in 
this sample are GJ 436b and HAT-P-20b, which have or^ 


bital eccen tricities of 0.15 and 0.016, respectively (Knut- 
son et al~1|2014b|. Despite this discrepancy, the data tor 


day-night circulation. 

We next consider whether or not the smaller planets in 
our sample show a different behavior. As noted earlier, 
HAT-P-19b, WASP-6b, WASP-39b, and WASP-67b are 
all best-matched by models with relatively efficient day- 
night circulation and are consistent with a range of at¬ 
mospheric metallicities. We identify two additional plan¬ 
ets in the literature with measured Spitzer secondary 
eclipse depths at 3.6 and 4.5 /im, masses greater than 
that of Jupiter, and predicted equilibrium temperatures 
less than 1200 K. As before, we compare the observations 
for these planets to our standard set of models, and plot 
the results in Figure [9| and flo) For WASP-80b (|Triaud 
et al.]|2015|, we find trie data require a model with rela¬ 
tively efficient day-night circulation and an atmosphere 
with at most a moderately enhanced metallicity. For GJ 
436b we obtain the best match using a model with effi¬ 
cient atmospheric circulation (this may be related to the 
planet’s non-zero orbital eccentricity) and a high atmo¬ 
spheric metallicity as discussed above. 

3.2.2. Temperature Inversions, Clouds, and Other Caveats 

Although the models described above suggest the pres¬ 
ence of some intriguing patterns, it is important to 
note that with measurements in just two bandpasses we 
cannot provide unique constraints on the atmospheric 
chemistries of these planets. In addition to the potential 
effects of disequilibrium chemistry and non-solar C/O ra¬ 
tios, the vertical thermal profiles of these planets can also 
alter the shapes of their dayside emission spectra. We ex¬ 
plore the potential effects of temperature invers ions us¬ 
ing a suite of models from Burrows et al. (2008), which 
allow us to artificially create atmospheric temperature 
inversions with varying strengths. Temperature inver¬ 
sions are created by adding a generalized gray absorber 
in the stratosphere above 0.03 bars with an absorption 
coefficient n of 0.1 cm 2 g -1 . A dimensionless parameter 
P n describes the efficiency of energy redistribution, with 
P n = 0.0 indicating recirculation on the dayside only, 
and P n = 0.5 indicating energy uniformly distributed 
between the two hemispheres. In these models the re¬ 
distribution is implemented via a heat sink located at 
pressures between 0.1 — 0.01 bars. 

As shown in Fig. B our measured eclipse depths for 
for four of the five planets presented in this paper require 
models with relatively efficient day-night circulation but 
we are unable to distinguish between models with and 
without a temperature inversion at a statistically signif¬ 
icant level. As with the previous set of models, WASP- 
lOb is best described by models with relatively inefficient 
day-night circulation. We note that the solar metallic¬ 
ity models with no inversion, which should be identical 
to the Fortney solar metallicity models, in fact display 
significantly stronger molecular absorption bands. This 
difference has been consistently observed in models Jor 
planets across a wide range of temperatures (e.g., |Dem-| 


ing et al. 120111 Tbdorov et al. |20 1 2 [ |20 1 3 [ |U ’ Itoui ke et al. 


these three planets suggest a clear pattern: planets with 
masses greater than that of Jupiter are best described 
by models with solar metallicity and relatively inefficient 


2014||Shporer et al.|2014|), and may be due to either the" 

use of different opacity tables, a change in their pressure- 
temperature profiles due to differences in their method 
for incorporating heat redistribution to the planet’s night 
side, or some combination of the two. 

We also consider the possibility that the dayside emis¬ 
sion spectra of these planets may be affected by the pres- 
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Fig. 11.— Solar metallicity atmosphere models showing the ef¬ 
fect of a temperature inversion Burrows et al.|2008|l. Model with 
efficient day-night circulation (P n = U.1J and either a. dayside tem¬ 
perature inversion (re = 0.1 cm 2 g —1 : light blue) or no inversion 
(re = 0.0 cm 2 g -1 ; purple) are shown as solid lines with band- 
integrated values over plotted as filled squares. All models assume 
an albedo of zero for the planet’s dayside. We also plot the corre¬ 
sponding models with inefficient day-night circulation ( P n = 0.4) 
and either a dayside temperature inversion (orange) or no inver¬ 
sion (red). Measured Spitzer secondary eclipse depths with their 
associated ler uncertainties are shown as black filled circles, with 
the exception of the 3.6 /im eclipse of WASP-67b where we plot 
the 2(j upper limit. 


ence of high altitude clouds or hazes, which have been ob¬ 
served in the transmission spectra of several cool gas gi- 


ant planets ( 

Line et al.12013 Jordan et al. 2013 

Knutson 

et al. 2014a 

INikoiov et al. 2015 

). A range of condensate 


(for a recent review see Marley et al. 2013), some of which 
may occur high enough up in the atmosphere to have 
a detectible effect. If the atmospheres of these planets 


have significant methane, then it may also be possible to 
form a layer of hydrocarbon soot via pho tochemistry in 
the upper layers of the atmosphere (Miller-Ricci Kemp-| 


ton et al. 2012 Morley et al. 


2013). Published optical 


transmission spectra lor WAbh’-bb indicate that it does 


indeed possess a high altit ude cloud layer (Jordan et al. 


2013 Nikolov et al. 2015), but there are no published 


transmission spectra available for HAT-P-19b, WASP- 
10b, WASP-39b, or WASP-67b. If these planets possess 
a cloud layer that is optically thick when viewed during 
secondary eclipse, we would expect such a cloud layer to 
suppress the strengths of the observed molecular absorp¬ 
tion or emission features, leading to a more blackbody¬ 
like spectrum. HAT-P-19b has a dayside emission spec¬ 
trum that is moderately inconsistent with that of a pure 
blackbody, which should produce a planet-star flux ratio 
that increases consistently towards longer wavelengths. 
WASP- 6 b, WASP-lOb, WASP-39b, and WASP-67b are 
closer to a blackbody shape, and we therefore cannot ex¬ 
clude the possibility that their dayside emission spectra 
are affected by clouds. However, we note that the transit 
technique is sensitive to more tenuous clouds and hazes 
than the secondary eclipse technique, as the slant op¬ 
tical path length traveled by light transmitted through 
the planet’s atmosphere is longer than the radial pat h 
traveled by light emitted by the planet (Fortney |2005). 
GJ 436b serves as a useful case in point, as it has an ap- 
parently featureless in frared transmission spectrum (e.g., 


Knutson et al. 2014a) but still exhibits strong spectra l 
features in its dayside emission spe ctrum (e.g., [Steven¬ 


son et al7]|2010 Lanotte et al.|[2014 ). 


3.3. Empirical Comparisons to the Sample of Observed 
Exoplanet and Brown Dwarf Emission Spectra 

Although our constraints on the atmospheric proper¬ 
ties of planets in individual cases are limited by the signal 
to noise of our measurements, it is possible that the ag¬ 
gregate sample may display trends that are not apparent 
on a case-by-case basis. We therefore examine the prop¬ 
erties of the planets in our sample using a series of model- 
independent metrics that allow us to construct a view of 
the global properties of this population. We begin by cal¬ 
culating the ratio of the 4.5 to 3.6 pm brightness temper¬ 
atures for each planet; this ratio is independent of both 
the planet-star radius ratio and the effective tempera¬ 
tures of the planet and star, and instead simply reflects 
the relative shape of each planet’s emission spectrum as 
determined by the presence of any molecular absorption 
or emission features. If we assume that molecular fea¬ 
tures are observed in absorption (i.e., none of the planets 
have dayside temperature inversions), we would expect 
the ratio of the 4.5 to 3.6 pm brightness temperatures to 
vary primarily as a function of atmospheric metallicity. 
In this scenario metal-rich planets will have relatively 
small brightness temperature ratios due to strong ab¬ 
sorption from CO at 4.5 pm and weak absorption from 
methane at 3.6 pm, while metal-poor planets will have 
correspondingly large ratios due to strong methane ab¬ 
sorption and weak CO absorption. 

We plot the measured brightness temperature ratios of 
the nine planets with published Spitzer secondary eclipse 
depths and effective temperatures cooler than 1200 K as 
a func tion of planet mass and equilibrium temperature 
in Fig. [T2| Interestingly, we see no evidence for a correla- 






























































TABLE 5 

Distance modulus, distance, apparent and absolute magnitudes for HAT-P-19b, WASP-6b, 
WASP-IOb, WASP-39B, and WASP-67b 
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Distance Modulus 
(m-M) 

Distance 

(pc) 

Apparent Magnitude 
m [3.6l m [4.5l 

Absolute Magnitude 
-^[3.61 -^4.51 

HAT-P-19b 

6.73 ± 0.25 

222 ± 28 

18.51 ± 0.23 

18.57 ±0.31 

11.79 ±0.53 

11.85 ±0.55 

WASP-6b 

6.65 ± 0.12 

214 ± 12 

17.84 ±0.22 

17.69 ±0.21 

11.19 ±0.38 

11.04 ±0.36 

WASP-10b 

5.58 ±0.22 

131 ± 14 

17.43 ±0.12 

17.09 ±0.12 

11.85 ±0.34 

11.50 ±0.33 

WASP-39b 

6.67 ±0.12 

216 ± 12 

17.80 ±0.19 

17.73 ±0.21 

11.12 ±0.33 

11.06 ±0.35 

WASP-67b 

6.64 ±0.14 

213 ±14 

> 18.32 

17.82 ±0.25 

> 11.68 

11.18 ±0.37 

600 700 800 
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Fig. 12.— Ratio of the 4.5 to 3.6 /mi brightness temperatures 
as a function of planet mass for the nine planets with predicted 
equilibrium temperatures cooler than 1200 K (colored points) or 
hotter than 1200 K (grey filled circles) and eclipse depths available 
in both Spitzer bands. Brightness temperature ratios are shown 
as filled circles (circular orbits) or squares (orbital eccentricities 
greater than 0.1; e = 0.30 for WASP-8b and e = 0.15 for GJ 436b) 
where the color corresponds to the predicted equilibrium tempera¬ 
ture of each planet at the time of secondary eclipse for the case of 
zero albedo and efficient redistribution of energy to the night side. 
Temperatures for the cool planets are t ake n from thi s study , and 
values for the hot planets are taken from|Hansen et al.|(|2014). The 
3.6 pm eclipse of WASP-67b and the 4.5 pm eclipse of GJ 4o6b are 
not detected at a statistically significant level, and we therefore plot 
the corresponding 2cr limits on the brightness temperature ratios 
for these planets. We note a general trend of increasing temper¬ 
ature ratio as planet mass increases; this trend appears to flatten 
out for planets more massive than Jupiter, consistent with the pre¬ 
diction that planets in this mass range should have envelopes that 
reflect the local metallicity of the protoplanetary disk. WASP-8b is 
an outlier in this plot, which may be related to its relatively large 
orbital eccentricity. 

tion with planet equilibrium temperature, but instead see 
tentative evidence for an increase in this ratio as planet 
mass increases from that of Neptune up to Jupiter. This 
is consistent with the hypothesis that lower mass plan¬ 
ets are more likely to have metal-rich atmospheres, while 
planets with masses comparable to or larger than that 
of Jupiter should have envelopes that reflect the (near¬ 
solar) metallicity of their local region of the protoplane¬ 
tary disk. We obtain a nearly identical plot if we simply 
take the ratio of the measured eclipse depths in these 
two bands, indicating that the varying effective temper¬ 
atures of the planets and their host stars have a relatively 
small effect on the measured eclipse depths as compared 
to changes in the observed molecular absorption or emis- 


Fig. 13.— Ratio of the effective dayside temperature derived 
from measured secondary eclipse depths to the irradiation temper¬ 
ature (predicted temperature at the substellar point for the case 
of zero albedo and no heat recirculation) as a function of planet 
mass for 51 transiting planets with published infrared secondary 
eclipse depths including the five new planets in this study. Planets 
with orbital eccentricities less than 0.1 are shown as filled circles 
where the color indicates their irradiation temperature, while plan¬ 
ets with eccentricities greater than 0.1 are shown as filled squares. 
The dayside and irradiation temperatures of these eccentric planets 
vary as a function of orbital phase, and is therefore not expected to 
follow a simple scaling relation. The predicted temperature ratio 
for the case of no redistribution of energy and zero albedo is shown 
as a dashed line, and the corresponding ratio for complete redis¬ 
tribution to the night side with zero albedo is shown as a dotted 
line. 


sion features. 

We evaluate the statistical significance of this trend 
by taking the error-weighted average brightness ratio for 
planets with masses less than 0.3 Mj up and comparing 
this ratio to the error-weighted average ratio for plan¬ 
ets with masses greater than 0.5 Mj up . We exclude GJ 
436b and WASP-67b from this calculation, as we only 
have lower and upper limits on their brightness ratios, 
respectively. We also exclude WASP-8 in light of its large 
orbital eccentricity, which may affect the observed prop¬ 
erties of its dayside atmosphere. We find that planets 
with masses less than 0.3 Mj up have an average bright¬ 
ness temperature ratio that is 2.1a smaller than the av¬ 
erage ratio for the more massive planets in our sample. 
We conclude more observations will be required in or¬ 
der to determine whether or not the observed trend is 
statistically significant. 

We next evaluate the relative efficiency of circulation 
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between the day and night sides of these planets by com¬ 
paring their observed versus predicted dayside emi ssion 
temperatures as des c ribed in Cowan & Agol (2011) and 


Schwartz & Cowan (2015). lor this exercise we con¬ 
sider" both the subset of planets with relatively low irra¬ 
diation levels including the five planets from this study, 
and also the full set of fifty transiting gas giant plan¬ 
ets with broadband infrared eclipse depths detected with 
greater than 3<r significance. We calculate the irradia¬ 
tion temperature, To = T* y/R*/a, for each planet based 
on published values for the stellar effective temperature 
and the planet’s scaled semi-major axis. This temper¬ 
ature corresponds to the predicted temperature at the 
planet’s substellar point assuming an albedo of zero and 
no redistribution of energy. The maximum hemisphere- 
averaged dayside effective temperature is then (2/3 ) 1 / 4 
of this value, and the temperature corresponding to com¬ 
plete redistribution of energy across the entire planet sur¬ 
face would be (1/4 ) 1 / 4 of this value. 

We calculate the effective dayside temperature mea¬ 
sured for each planet using all available infrared (A > 
1 /im) eclipse depth measurements. We convert the set 
of individual dayside brightness temperatures across all 
bandpasses into a combined dayside effective tempera¬ 
ture using two different schemes: error-weighted mean 
(wavelengths with small uncertainties contribute most 
to the estimate) and power-weighted mean (wavelengths 
that capture large fractions of the planetary emission 
contribute most to the estimate). The uncertainty in 
the averaged dayside effective temperature and in its 
ratio with the planet’s irradiation temperature is then 
estimated using a 1000-step Monte Carlo method with 
equa l numbers of trials drawn from each scheme (see 


Schwartz & Cowan 2015 for more details). We use the 
ratio of the measured vs predicted (irradiation) dayside 
temperatures to determine whether or not the correla¬ 
tion between planet mass and day-niglrt circulation effi¬ 
ciency inferred in our cool planet sample is consistent 
with the larger sample of published secondary eclipse 
data for transiting gas giant planets with a range of in¬ 
cident flux levels. We plot this ratio as a function of 
planet mass in Fig. |13[ and find that our new observa¬ 
tions clearly indicateuie existence of a correlation be¬ 
tween planet mass and day-night circulation efficiency 
for planets with lower irradiation levels. This correlation 
also appears to be consistent with the published measure¬ 
ments for more highly irradiated planets, although there 
are fewer low-mass planets in this sample. We conclude 
that lower-mass planets may be more reflective and/or 
better at transporting heat to their night sides. Although 
the current data do not allow us to distinguish between 
these two hypotheses, the combination of visible-light 
secondary eclipse observations and infrared phase curve 
measurements could provide a definitive answer to this 
question in future studies. 

Lastly, we create a color magnitude plot for the five 
planets in this study that allows us to compare them to 
ultracool brown dwarfs and directly imaged planets span¬ 
ning the same range of effectiv e tem peratu res using the 


same procedure a s described in Triaud (2014a) and Tri- 


aud et al. (2014b). In order to maintain consistency with 


the other planets in this diagram, we obtain stellar pa¬ 
rameters for the five planets presented in this paper from 



Fig. 14.— Diagram comparing the 4.5 fim absolute magni¬ 
tudes and 3.6—4.5 ^m colors of HAT-P-19b, WASP-6b, WASP- 
10b, WASP-39b, and WASP-67b (blue filled circles and square) to 
those of previously published planets a nd brown dwarfs in th e same 
magnitude range (figure adapted from Triaud et al. |2014b[ ). Pub- 
fished planets are plotted as grey circles, while ult ra-cool dwarfs 
with s pectral types between M5—Y1 obtained from|Dupuy & Liu| 
(2012|) are plotted as diamonds with colors indicating their spectral 
type (orange corresponds to hotter stars and brown indicates cooler 
stars). WASP-67b is not detected in the 3.6 ^m band and is there¬ 
fore is represented with a blue square and arrow to indicate the 
corresponding 2 a lower limit on its 3.6—4.5 f-im color. The black 
line indicates the color of a 1 Rj up object emitting as a blackbody 
over a range of effective temperatures less than 4000 K. 


TEPCAT (Southworth 2011) and c ompute their abso¬ 
lute magnitudes u s ing re lations by Flower (1996) and 


Torres & Gimenez (2010). Distance moduli were calcu- 
lated by comparing them to visual magn itudes, which we 


take from the UCAC4/APASS catalog (Zacharias et al. 


20131. Distance moduli, derived distance, visual and ab¬ 


solute magnitudes are provided in Table 5] and visually 
represented in Figure [14] We calculate these qua ntities 


using WISE 1 and WISE 2 apparent magnitudes (Cutri 


et al.||2013l, which a re a good match for t he Spitzer 3.6 
and 4.5 gm bands ( |Triaud et al. 2014b), and list the 
resulting values in Tablets} 

As shown in Fig. [HJ allfive planets have similar mag¬ 
nitudes in this diagram, albeit with varying colors. A long 


with GJ 436b (Lanotte et al. 12014 1 and WASP-80b (Tri- 
aud et aI7||2015 ), these five planets have the smallest ab¬ 


solute magnitudes of any extrasolar planets observed to 
date and demonstrate the benefit of extending the cur¬ 
rent Spitzer observations to smaller and cooler planets. 
Three of the five planets appear to be bluer than a black- 
body with the same absolute magnitude, placing them 
near the late L dwarfs in this diagram. WASP- 6 b and 
WASP-67b lie on or near the blackbody line, where the 
2cr upper limit on the 3.6 pm eclipse depth for WASP- 
67b indicates that it is al most certain ly redder than a 
blackbody. As discussed inlTriaud et al. (|2014b|), irradi¬ 


ated exoplanets have optical properties that appear to be 
more diverse than those of self-luminous brown dwarfs, 
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as demonstrated by the varying colors of the planets in 
our sample. Although part of this variation likely reflects 
intrin sic variations in surface gravities and met allic ities 
(e.g., Zahnle & Marley 2014), our analysis in 13.2 sug¬ 
gests that disequilibrium chemistry, cloud formation, and 
variations in atmospheric circulation patterns may also 
contribute to the observed diversity of dayside emission 
spectra for the transiting planet sample. 

4. CONCLUSIONS 

We present new secondary eclipse depth measurements 
for HAT-P-19b, WASP-6b, WASP-lOb, WASP-39b, and 
WASP-67b in the 3.6 and 4.5 fj, m Spitzer bands. We 
use our measured center of eclipse times to constrain the 
orbital eccentricities of these planets and find that HAT- 
P-19b, WASP-6b, WASP-39b, and WASP-67b all appear 
to have circular orbits. Our new secondary eclipse times 
confirm the previously reported non-zero orbital eccen¬ 
tricity for WASP-lOb, and we carry out a joint fit to all 
of the available transit, secondary eclipse, and radial ve¬ 
locity data in order to provide improved constraints on 
this planet’s orbit. 

We combine our observations with published secondary 
eclipse depths for other transiting gas giant planets and 
find that the most massive planets are best matched by 
models with low albedos and relatively inefficient day- 
night circulation, while smaller planets are best matched 
by models with higher albedos and/or efficient day-night 
circulation. Previous analyses of the availabl e secondary 
eclipse data for transiting hot Jupiters (Cowan fe Agol 


2011 |Perez-Becker &; Showman 2013; Schwartz & Cowan 


20151 have found that planets with lower irradiation lev¬ 

els appear on average to have more efficient recirculation 
of energy than their more highly irradiated counterparts, 
but did not note any correlation with planet mass. Our 
results are consistent with our conclusions in 1 Schwartz! 


& Cowan (2015) and Wong et al. (20151, where we ex- 
amined the sample of planets with infrared phase curve 
observations and found that the two most massive plan¬ 
ets in this sample appear to have systematically lower 
albedos than their less massive counterparts as inferred 
from a simple energy balance calculation. In the future, 
measurements of the secondary eclipses of these planets 
at optical wavelengths will help to establish whether or 
not the observed variations are primarily due to changes 
in the Bond albedos of these planets or to differences in 
their atmospheric circulation. 

All of the planets in this study have predicted equi¬ 
librium temperatures cooler than 1200 K, and we there¬ 
fore expect that the ratio of methane to CO in their 
atmospheres should vary as a function of their atmo¬ 
spheric metallicity. When we combine our observations 


with published secondary eclipse data for four additional 
planets in this same temperature range, we find that 
planets with masses greater than that of Jupiter are 
best matched by relatively low metallicity atmosphere 
models. Although our results are less well constrained 
for the smaller planets in our sample, we examine their 
aggregate properties by taking the ratio of the 3.6 to 
4.5 fim brightness temperatures and plotting this ratio 
as a function of planet mass and temperature. We find 
that there is no detectible correlation between this ratio 
and the equilibrium temperatures of these planets, but 
we do see tentative evidence for a trend in mass consis¬ 
tent with the existence of an inverse correlation between 
planet mass and atmospheric metallicity. We note that 
with just two wavelengths of data our ability to constrain 
the atmospheric metallicities of these planets is degen¬ 
erate with other atmospheric properties, including the 
presence or absence of dayside temperature inversions as 
well as the potential effects of high altitude cloud layers. 

Future observations of the dayside emission spectra of 
these planets spanning a broader range of wavelengths 
will be able to determine whether or not the observed 
trends are in fact due to variations in atmospheric metal¬ 
licity as suggested here. Although these planets are too 
cold and faint to have detectible secondary eclipses in the 
near-infrared wavelengths accessible to the Hubble Space 
Telescope and ground-based telescopes, future spectro¬ 
scopic observations at longer wavelengths with the James 
Webb Space Telescope will provide invaluable informa¬ 
tion about their atmospheric compositions. For planets 
without high altitude clouds or hazes, measurements of 
water absorption in transit with WFC3 on HST can also 
provide independent co nstraints on the inferred atmo¬ 
spher ic metallicity (e.g., Kreidberg et al.||2014 Benneke 
20151. In the short term/observations of additional plan¬ 


ets in the Spitzer bands will help to determine whether 
or not the observed correlation between planet mass and 
3.6 to 4.5 /im brightness temperature ratio is statisti¬ 
cally significant, and will allow us to better understand 
the differences between cool transiting gas giant planets 
and the sample of directly images planets and ultra-cool 
brown dwarfs with similar temperatures. 
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